Electron heating in metallic resistors at sub-Kelvin temperature 
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In the presence of Joule heating, the electronic temperature in a metallic resistor placed at sub- 
Kelvin temperatures can significantly exceed the phonon temperature. Electron cooling proceeds 
mainly through two processes: electronic diffusion to and from the connecting wires and electron- 
phonon coupling. The goal of this paper is to present a general solution of the problem, in a form 
that can easily be used in practical situations. As an application, we compute two quantities that 
depend on the electronic temperature profile: the second and the third cumulant of the current 
noise at zero frequency, as a function of the voltage across the resistor. We also consider time 
dependent heating, an issue relevant for experiments in which current pulses are used, for instance 
in time-resolved calorimetry experiments. 
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PACS numbers: 72.15.Lh,73.50.-h,73.23.-b,73.50.-(-d,72.70.+m 



I. MOTIVATIONS AND OUTLINE 

When performing electrical measurements, the signal 
to noise ratio can usually be improved by simply increas- 
ing the currents or voltages. In low-temperature exper- 
iments, this procedure is problematic because of Joule 
heating, which can affect the temperature of the cir- 
cuit under investigation or the temperature of resistors 
on bias lines, leading to excess noise. Particularly criti- 
cal is the situation in the sub-Kelvin range, because the 
temperature of the electrons decouples from the lattice 
temperature^"^. A very conservative, wide-spread rule 
of thumb among experimentalists is that the voltage V 
across the small conductors should not exceed fesTph/e, 
with Tph the lattice (phonon) temperature. In contrast, 
macroscopic components, like commercial resistors, are 
believed to be immune to electron heating. In fact, the 
first rule is severe, and the second assertion is often in- 
correct. The goal of the present article is to provide the 
experimental physicist with easy evaluation tools of heat- 
ing effects, in order to optimize experiments. 

The important parameters are the voltage V, the re- 
sistance R, the lattice temperature Tph, the resistor vol- 
ume J7 and a parameter S that describes electron-phonon 
coupling. The first step is to calculate the characteristic 
temperature Ts, which is the temperature that the elec- 
trons would reach if cooling would occur only through 
the coupling to a bath of zero-temperature phonons: 
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(1) 



The average electron temperature can then be directly 
read from Fig. 2, the central result of this work, in which 
the voltage V, the average temperature Tav and the lat- 
tice temperature Tph are all expressed in units of Ts . In 
the section IIA, we explain how this result is obtained, 



and give analytical expressions in various limits. The 
results are used to calculate the second and third cumu- 
lant of the current noise produced by the resistor (section 
IIB, Fig. 4). Numerical applications are carried out ex- 
plicitely in section IIC, showing in particular that heat- 
ing in commercial resistors can be important. In section 
III, we address time-dependent situations, and calculate 
how fast electrons heat up in a resistor when a current is 
applied, and how fast they cool down when the current 
is switched off. For small voltages {eV ^ ksT^), the 
variations of temperature in both transients is exponen- 
tial, with the diffusion time td across the whole conduc- 
tor as a characteristic time (see Fig. 5). In the opposite 
limit {eV ^ ksTs), heating is exponential, but cooling 
proceeds very slowly, with a powerlaw dependence (see 
Fig. 6). The timescale is the electron-phonon scattering 
time Tc-ph at temperature Ts, defined by Eq. (17). As 
a numerical application, we consider in section IIIC a 
situation where repeated current pulses are applied to a 
resistor, and compute the time-dependence of the elec- 
tron temperature (Fig. 7). 



II. STATIONARY SITUATIONS 



A. Solution of the heat equation 



When a voltage V is applied to a two-terminal resis- 
tor (see top of Fig. 1), the Joule power V'^/R is deliv- 
ered to the electrons. This power can dissipate by two 
mechanisms: the first one, which dominates at room tem- 
perature or in macroscopic resistors, is phonon emission. 
It follows, at temperatures well below the Debye tem- 
perature, a T"(a;) — Tpj^ dependence, with 4 < n < 6, 
and T{x) the local electron temperature, Tph the phonon 
temperature^. The second mechanism is the simple diffu- 
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sion of the energetic electrons out of the resistor. The en- 
ergy is then dissipated in the connecting leads, which are, 
in typical situations, large and low-resistive. The balance 
between the Joule power and the two cooling mechanisms 
can be expressed in the form of a heat equation^ 

with X the position along the resistor in reduced units 
(x runs from to 1), Lo = TT^k%/3e'^ the Lorenz num- 
ber, the resistor volume, S the clectron-phonon cou- 
pling constant (typically S ~ 2 nW//xm'^/K^ for good 
metals''). The left hand side of Eq. (2) accounts for heat 
transport by electron diffusion, which is expressed by the 
Wiedemann-Franz law, stating that the electron thermal 
conductivity is proportional to the product of the elec- 
trical conductivity and the electron temperature. The 
following assumptions have been made to write Eq. (2): 

1. The electron temperature T{x) is assumed to be 
well defined locally, i.e. the local electron energy 
distribution function is a Fermi function. This re- 
quires that the thermalization of electrons among 
themselves (e.g. by Coulomb interaction) occurs 
faster than the diffusion of electrons across the 
resistor^, a condition usually obeyed except for 
short wires (length < 50 fim) made of very pure 
materials^. 

2. The last term of the equation, which describes cool- 
ing by phonons, assumes that the lattice temper- 
ature Tph does not depend on the local electron 
temperature T{x). Corrections due to the Kapitza 
resistance between the phonons of the resistive film 
and the substrate could in principle be included^'^°, 
but their contribution is not essential in practice. 

3. The heat power transfered to phonons was taken 
proportional to T"(.t) — T^^ with n = 5. Theoret- 
ically, it is predicted that the exponent n, which 
is related to a £^^~" dependence of the electron- 
phonon scattering rate with electron energy E, 
can range from 4 to 6, depending on the rela- 
tive sizes of the thermal phonon wavelength and 
the electron mean free path, on the dimension- 
ality of the phonon system, and on the dynam- 
ics of impurities^. In most experiments, values 
close to n = 5 have been found (see discussions in 
Refs. 5,7,11), therefore our choice. The calculations 
can however be easily extended to other values of 
n, and the results presented in Fig. 2 apply, with 
another definition of Ts, as discussed in the follow- 
ing. 

4. Radiative cooling^ ^, which has a negligible effect in 
resistors connected to large, non-superconducting 
contacts^^, is neglected. 



The heat equation (2) has to be solved with boundary 
conditions for T{x) at x = and x = 1. When the con- 
necting wires to the resistor are low-resistive and very 
large compared to the resistor, as is the case for macro- 
scopic resistors made of thin and narrow metallic stripes 
of metal, one can assume T(0) = T{1) = Tph- This sim- 
ple hypothesis will be made in the following. For on-chip 
thin-film resistors, heating of the contact pads themselves 
may however not be negligible^"^. 

Before a general solution of Eq. (2) is presented, we re- 
call simple limits. The so-called interacting hot-electron 
limit^ is obtained by neglecting phonon cooling: 

T(.) = \/tp\+^.(1-.)(|^)' (3) 

(see left panel of Fig. 1, dashed lines). For Tp^ = 0, 
the maximal temperature is r(|) = {V^/2n) eV/ks — 
0.28eV/kB, and the average temperature is Tav = 
J^T{x)dx = {VS/8)eV/kB ^ 0.22eV/kB. Electron- 
phonon coupling further reduces the temperature, so that 
this is an upper bound on the average electron tempera- 
ture, which numerically reads, keeping^^ now Tph: 

In particular, the rule of thumb eV = ksTph corresponds 
to a 2.5% average overheating of the electrons. Numeri- 
cally, one obtains the equivalent expression 

Tav [mK] < ^Tp\[mK]+(2.5xl/[MV])^ (5) 

In the opposite limit where cooling by diffusion can be 
neglected, the electron temperature is homogeneous and 
equal to^ 

T=(Tp\ + T|)'/^ (6) 

In the following, we call this limit "the fully thermalized 
regime" . In the limit Tph = 0, the temperature grows as 

(sec Fig. 3). 

In intermediate regimes, the temperature profile is ob- 
tained by solving numerically the heat equation. For gen- 
erality, it is convenient to rewrite it in reduced units. 
One possibility is to take as a reference the "cross-over 

temperature"* Tco = (T,QRe'^ / kg) which is the en- 
ergy scale for which phonon cooling and diffusion cooling 
are equaly important. This temperature is an intrinsic 
quantity for the resistor, which in particular does not de- 
pend on V or Tph- However, it does not correspond to 
the electron temperature in any limit, therefore we pre- 
fer to take as a reference Ts, keeping in mind that it 
depends on V. Defining 9{x) = T{x)/T^, v = eV/ksT^ 
(= (ey//cBTeo)3/5), V = Tpi,/T^ (oc F-2/5), Eq. (2) 
reads: 

^e^{x) = ^v^9'{x)-el,-i). (7) 
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The temperature profile being symmetric with respect 
to the middle of the wire, Eq. (7) needs to be solved for 
< X < ^ with the boundary conditions 9{0) — 9ph and 
6''(i) = 0. Instead of solving this non-linear differential 
equation with boundary conditions specified at different 
points, it is convenient to rewrite it in the following in- 
tegral form: 



— n ^ = / 

(8) 

with 9^ = 6'(i) and A = [91^ + 1 j . The value of 9m 

is obtained by solving Eq. (8) for x — 1/2. In the left 
panel of Fig. 1, the temperature profile along the resistor 
is given for i; = 1, 2, 3, 10 and 30, assuming Tph — 0. At 
w < 1, one recovers the result of Eq. (3), plotted as dashed 
lines: phonon cooling can be neglected. For v > 10, it 
is an excellent approximation to take 9m = A in Eq. (8), 
a value that does not depend on v. T(x) is then a func- 
tion of XV only, which is essentially constant around the 
middle of the wire, whereas at a distance b/v from the 
contacts one obtains the profiles shown in the right panel 
of Fig. 1 for 6'ph — (see note 15), 0.5 and 1. The charac- 
teristic length over which the electron temperature varies 

from Tph to Ts is therefore is = L/v = ^i^/^if^ph, 
where the "electron-phonon length" 
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which V is changed at a fixed Tph, since 0ph oc V~'^/^. 
In order to visualize how temperature increases with V 
at a given phonon temperature, we plot in Fig. 3 the av- 
erage temperature Tav(T^) for various Tph, with V and 

Tph given in units of* T^o = (srii?e^/fc^) , which 
is constant for a given resistor. The range in voltage 
V is the same as in Fig. 2. We have used the relations 

Z;5/3 = eV/kT,,- (Tav/Ts) 1-2/3 = Tav/T^. 



ph 
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{p being the resistivity) is defined, foUowing-'^^ Ref. 6, as 
the resistor length for which AkBTs/R is equal to the 
current noise in the interacting hot-electron regime -^el. 
This length is typically of the order of a few fim for volt- 
ages of the order ot 1 mV. 

From the complete temperature profile 9{x), the av- 
erage electron temperature Tav is obtained using Tav = 
/q 9{x) dx. The central result of this work is the re- 
sulting plot, shown in Fig. 2, of the average temperature 
Tav as a function of the voltage V, both in units of Ts, 
for Tph/Ts = 0, 0.25, 0.5, 0.75 and 1. At 6'ph = and 
I' < 2, Tav ~ (dotted line), whereas for u > 4, 

7av/Ts w 1 — 1.16/i; (dashed fine). This l/v dependence 
is due to the crossover regions of width ~ 5/v at the re- 
sistor ends. Figure 2 can be directly used to read out the 
average electron temperature Tav for a given set of exper- 
imental parameters (V, Tph), after having computed Ts 
with Eq. (1). Interestingly, the corresponding curves for 
other exponents of the temperature in the last term of 
Eq. (2) (n = 4 or 6 instead of 5) are almost identical (see 
inset), and the same curves can be used to evaluate Tav, 
however with the generalized definition of the reference 

temperature Ts = {V^/^nR}^^'' . 

However, because of the use of reduced units which de- 
pend on V, the w-dependence of Tav at a fixed value of 
9ph shown in Fig. 2 does not correspond to a situation in 



B. Second and third cumulants of the current noise 
at zero frequency 

The temperature profile can be used to evaluate the 
current noise properties of the resistor. We focus here 
on the second (S'2) and third {S3) cumulants of noise at 
low frequencies {Huj <C eV for S2, fu^ eV, h/T£) for S3 
(see Ref. 17 ), with td — /D the diffusion time, D the 
diffusion constant): 

^2 = ^JJ'^t ('5/(0)«(i)) 



and 



^3 = dtidia («(0)<5/(ti)«(i2)) 



1 


resistance R, volume n 


1 

ph 


T(x) 




FIG. 1: (Color online) Top: Resistor biased by a voltage V 
and placed between two connecting wires in which the elec- 
tron temperature T and the phonon temperature Tph are 
equal. Left panel: Solid lines: temperature profile in the 
resistor for different values of v = eV/ksTY^ with Te = 
(l/^/Efii?) Dashed lines: temperature profile expected 
when phonon cooling is neglected (Eq. (3)). Right panel: tem- 
perature profile near the ends of the resistor for n > 10 and 
6'ph = Tph/TE = 0, 0.5, 1. 
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with 51{t) = I{t) — (/) . It has been shown that when 
phonon coohng can be disregarded and eV ^ k^Tpi^, S2 
and S3 are proportional to the appUed current^' 5*2 = 
F2 X 2eJ and S3 = F3 x e^I with F2 and F3 generahzed 
"Fano factors" . When furthermore the rate of electron- 
electron interaction is negligible compared to 1/td, the 
distribution function is not a Fermi function, but a 
function with two steps^'^^, and^^"^° F2 = ^ « 0.33; 
F3 — Ri 0.067. In the opposite limit, where electron- 
electron interaction is strong, electrons thermalize locally 
to distribute in energy according a Fermi function, and 
the temperature profile is given by Eq. (3). One then 
obtains22.35 = ^ « 0.43 and^^ i^g = - A « 0.248. 

In presence of strong phonon cooling {v ^ 1), the elec- 
tron temperature becomes homogeneous, at a value 
smaller than eV/ks- It is then expected^ ^ that F2, F3 — + 
0. 

For intermediate coupling to phonons, F2 and F3 de- 
pend on the voltage across the resistor. Their value 
are obtained from the full solution of the heat equation: 
the second cumulant is given by a Johnson-Nyquist-like 
formula^ ^2 = 4fcBTav/i? in which the noise tempera- 
ture is the average electron temperature Tav, yielding 
F2 = 2kBTa^/eV = 26'av/w. This formula can be un- 
derstood as resulting from the added Johnson-Nyquist 
noise of small sections of the resistor, each at a temper- 
ature T{x). The decay of F2 at large V was discussed in 
Ref. 22 , and the complete crossover was calculated in 
Rcf. 23 by numerical integration of Eq. (2). In turn, F3 




V = eV / k T 



FIG. 2: (Color online) Average temperature Tav in units of 
Te, as a function of 1; = eV/ksTs, for 9ph = Tph/Ts — 0, 
0.25, 0.5, 0.75, 1 (from bottom to top). The dotted line corre- 
sponds to the low-u approximation Tav = 1 the dashed- 
dotted line to the large-u approximation Tav ~ 1 — 1.16/t;. 
Inset: at dph = Tph/Ts = 0, comparison of the evolution of 
the average temperature with v for various exponents n of 
the temperature in the expression of the heat flow through 
electron-phonon coupUng. 



is given by (see Appendix) 

36 1 

F3 = — dxdy—Gi{e,x,y){e{y)-2e,^} (10) 
Jo '^\^) 

where Gi{6,x,y) the Green's function such that 
iy^ + ^vH\x))Gi{e,x,y) = 6{x - y) and 
Gi(6',0,y) = Gi{e,x,Q) = 0. The calculation of F3 
is detailed in the Appendix. The right panel of Fig. 4 
shows the voltage dependence of F2 (blue line) and 
F3 (red line) as a function of v (bottom axis) and 

5 /2 

L/Lo-ph — (\/3w/8) (top axis), for Tph = 0. Also 
shown with a dashed line is the curve obtained for F2 
when electron diffusion is neglected^^, using Eq. (6) 
(dashed line), which gives F2 = 2/v (x V~^^^. In turn, 
at large voltages, F3 (x (x If one considers 

a situation where the resistor length L is varied at 
constant current, then F2 cx 1/L and F3 cx 
The decay to zero of F2 and F3 from the interacting 
hot-electron values (-^ and ^ — ^) is therefore very 
slow, as already pointed in Ref. 21,23. The non-Gaussian 
character of the current noise, evidenced by F3 ^ 0, is 
washed out at L/L^-ph ^ 10. 

C. Examples 

We consider here a few cases illustrating the use of 
the results given in the preceding sections. As a first 
example, we consider a 10 //m-long, lOOnm-wide and 
5 nm-thick Cr resistor with resistance R = 25 kfl like 
these used in Ref. 25, biased at F = 1 mV and placed at 




eV / k T 

B CO 



FIG. 3: (Color online) Average temperature Tav as a func- 
tion of voltage V, for various temperatures Tph, all in units of 

Tco = {T,D,Re'^/kiy^^'' . The value of Tph/Tco is given by the 
intersection of the curves with the vertical axis. Blue dashed 
line is the reference temperature Te. 
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Tph = 20 mK. Assuming Scr = 2nW/^ni3/K5, the char- 
acteristic temperature is = 1.3 K, and the voltage and 
phonon temperature in reduced units v = eV/kBT^ ^ 9 
and 0ph = 0.015. The noise temperature is directly read 
from Fig. 2: Tav 0.87 x ^ 1.15 K. At this vohage, 
heating of the resistor is thus very important, an effect 
which hindered the authors of Ref. 25 from drawing clear- 
cut conclusions from Coulomb blockade measurements 
at finite voltage. Increasing the resistor volume fl with 
"cooling fins" can help decreasing electron heating^, but 
such a procedure is extremely inefficient since the char- 
acteristic temperature Ts decreases as fl^^^^ only. 

As a second application, we now consider a com- 
mercial macroscopic surface mount resistor, with R = 
500 n. Such resistors, made of thin (~ 10 nm) NiCr 
films^^ with resistivity p ~ 100 /if2 cm and dimensions 
'--^l mmxO.2 mm, were used as bias resistors in mea- 
surements of the state of superconducting Josephson 
Q-bits^^'^* performed at 15 mK, with a bias current 
~ 0.8 /iA, resulting in a voltage V ^ 400 //V. The 
corresponding temperature scale Ts ~ 150 mK yields 
V — eV/kBTj: ^ 30 and 0ph ^ 0.1, hence, from Fig. 2, 
Tav ~ ~ 150 mK. Even in such a macroscopic resis- 
tor, the volume is not sufficient to provide with enough 
electron-phonon coupling, and heating is important. In 
the next section, we show how this heating is limited 
when pulses are used instead of static voltages. 




V = eV / k T 
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FIG. 4: (Color online) Solid lines: Fano factors F2 = 52/26/ 
(blue) and F3 = S's/e^l (red) for the zero-frequency second 
and third cumulant of noise, from the hot electron limit to 
the fully thermalized regime. Dashed-dotted line: asymptotic 
dependence of F2 neglecting electron diffusion: F2 — 2/v. Top 
axis is resistor length over Le-ph (see Eq. (9)). 



III. TIME-DEPENDENT SITUATIONS: 
SWITCHING ON AND OFF JOULE HEATING 



The case of a constant voltage V across the resistor, 
which was investigated above, can be extended to the 
case of slowly varying voltages directly. However, when 
V changes on timescales shorter than the diffusion time 
or than the electron-phonon scattering time^^ (see be- 
low), the previous results cannot be used to calculate 
instantaneous temperatures. These issues are solved by 
adding to the heat equation (2) a time-dependent term 
dQ/dt = CedT/dt, with Cg = jT^l the electronic heat ca- 
pacity, 7 = {tt^ /3)kgVp (from Fermi liquid theory), with 
h'p the density of states at Fermi energy (spin degeneracy 
included). When V{t) = Vf{t), the time-dependent heat 
equation can then be rewritten, in reduced units, as 
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fir)) (11) 



where r = t/ro is the reduced time with td = L^/D 
the diffusion time. Note that the reference temperature 
Ts used to define 9{x) ~ T{x)/Ty: is calculated with 
the voltage scale V, not with the time-dependent value 
V{t). In the following, we treat more explicitely two sit- 
uations: how a resistor heats up when the voltage is ap- 
plied at t = 0, i.e. /(r) = H{t), and how a resitor 
cools down when the voltage is set to zero at t = 0, i.e. 
/(r) = 1 — H{t). Here, H{t) is the Heaviside function 
(0 for r < 0, 1 for r > 0). These situations also allow to 
describe experiments in which current or voltage pulses 
are used like, for example, when measuring the switching 
rate of Josephson junctions^^. Understanding how the 
pulse characteristics can reduce the noise in such mea- 
surements is therefore important to design the readout 
of superconducting Q-bits. 

When a voltage is applied, the linear drop of the elec- 
trical potential, which results from the collective charge 
modes, establishes after an RC time, where the capaci- 
tance C is the capacitance of the wire to ground. This 
time is generally much shorter than the time necessary 
to build up the temperature profile, which involves dif- 
fusion of individual electrons. Hence, we consider here 
that Joule heating is homogeneous as soon as a voltage 
is applied. When v < 1, the temperature profile is en- 
tirely determined by the temperature at the ends of the 
resistor, therefore a steady-state regime is reached only 
when the electrons have diffused across the whole resistor 
and the characteristic time is the diffusion time td. If 
w ^ 1, the transient is shorter because, apart from very 
close to the ends, the temperature is mostly determined 
by a local equilibrium between Joule heating and phonon 
emission. We now treat quantitatively these two limits. 
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Small V limit 



If Tj < 1, 0{x,t) — 9p\i <C 1 even when the stationary 
regime is reached, and Eq. (11) reduces to 



.2n2 



dr dx^ 



(12) 



As shown in section IIA, the proper energy scale when 
?; < 1 is eV, and the solution of Eq. (12) that satisfies the 
boundary conditions T(0, r) — T{1, r) = Tph reads 



kBT{x, t) 



eV J \ eV 
with afc(T) solution of 



's — ph ^ I afe(r) sin (vrfcx) 

(13) 



k odd 



dafc(r) 
dr 



(14) 



In the case where /(r) = H{t), heating is then given 
by afc(T) = (24/7r^P) (1 - e""^''^'^), and in the case 
/(r) 1 — H{t), cooling from the profile (3) follows 
o.k{T) ~ (24/7r^fc^) e^'^ ^ . Corresponding temperature 
profiles at various times are plotted in the top panels 
of Fig. 5 assuming Tph = 0, whereas the time evolution 
of the average temperature Tav is plotted in the bottom 
panels. At very short times, the average temperature 
grows as ksTav/eV = \/6t/tt. At t > 0.01, a better 
approximation is 



V3 



eV 



Vl-exp(-lOr), 



(15) 



which cannot be distinguished from the exact solution in 
Fig. 5. At T ~ 0.5, the asymptotical temperature profile 
given by Eq. (2) is essentially established. Conversely, 
after the voltage is turned off, the temperature decay is 
well approximated (within the line width in Fig. 5) by 



V3 



eV 



exp(— 5 r). 



(16) 



Hence, in the low-voltage regime, heating and cooling 
occur exponentially, and the timescale is the diffusion 
time td ■ In the example of the 500 commercial resistor 
in section IIB, tjj ~ 3 ms. For metallic wires made of pure 
materials^^ with an elastic mean free path of the order 
of 40cm^/s, Td ^ 20 ns for a length L ^ 20 /im. 



B. Large v limit 

If ti ^ 1, it was shown in section IIA that the tem- 
perature becomes almost homogeneous in the wire. The 
relevant timescale is then the electron-phonon scattering 
time""^^'^"''^^ at the characteristic temperature Ts : 



To-ph (Je) 
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3 



D 



(17) 



with Ls the characteristic length for the variation of 
T{x) introduced in section IIA. Numerically, 7/E « 
0.03 /isK^. Using = 2t* with 



r* = t/Te-ph (Ts) , 



Eq. (11) reduces to 
902 



dr* 



which for 0ph = / = is simply equivalent to 

dT _ T 
'dt " "re_ph (T) 



(18) 



(19) 



(20) 



expressing that the instantaneous decay rate of T is ex- 
ponential with a characteristic time Tc_ph (T) . Equa- 
tion (19) yields 



dw 



ph 



,,5/2 



2r* 



(21) 



When /(r) = 1 — H{t) and Tph — 0, the temperature 
decay from T-^ has a simple form: 



(1 -|-3r* 



-1/3 



(22) 



This temperature decay, which was directly measured 
in Ref. 11 , follows a power law only, so that it takes 




FIG. 5: (Color online) Time evolution of the temperature pro- 
file (top panels) and of the average temperature Tav (bottom 
panels) in the limit v <^ 1, for Tph = 0. Left panels: heat- 
ing sequence of the resistor when V{t)/V = H{t) (as shown 
in inset); right panels: cooling sequence of the resistor when 
V{t)/V = 1 — H{t) (as shown in inset). The profiles are 
plotted at various values of t/ro, with td ~ L^/D the dif- 
fusion time. The colors of the solid curves in the top panels 
corresponds to those of the square symbols in the bottom pan- 
els. The curves in the bottom panels cannot be distinguished 
from ksTav/eV ~ ^1 " exp(-lOT) and ^ exp(-5r), re- 
spectively. 
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a very long time to recover the base temperature after 
the voltage is set to 0, which is due to the divergence 
of To_ph (T) when T 0. The results of Eq. (21) with 
/ = 1 (heating) and f — (cooling) are plotted in Fig. 6 
in the case Tph — 0, with linear (top) and logarithmic 
(bottom) time scales. The temperature rise is well ap- 
proximated by 0{t*) « y/2 T* when t* < 0.2 (dashed 
line) and 9{t*) « 1 - 0.86 exp(-4.2 r*) when r* > 0.2 
(dotted line). More generally, when r* <C 1, for Tph 7^ 0, 

Oir*) « (0) + 2(1 + 0^^- 05 (0))r*. (23) 

Even though w 3> 1, we now estimate cooling by elec- 
tron diffusion to the connecting leads. Starting from 
a constant temperature Tq, cooling by diffusion follows 

Eq. (13) with atir) - (^)' and T^v/To - 

exp(— 5r). Because of this exponential dependence, to 
be compared with the powerlaw (22), diffusion can con- 
tribute to the cooling when t becomes comparable to tjj . 



C. Numerical application 

In experiments where the voltage is applied in repeated 
pulses, heating is reduced, and the temperature oscil- 
lates in time. To illustrate this effect, we reconsider 
the second example of section IIC, but we now assume 
that the voltage is applied during short pulses of du- 
ration tp = 0.1 /is, repeated every period t,, = 20 /is 




t / T (T ) t / T ,(T ) 



FIG. 6: (Color online) Evolution of the temperature with time 
in linear (top) and log (bottom) scale in the limit « 2> 1, for 
Tph = 0. Left panels: heating sequence of the resistor when 
V{t)/V = H{t) (as shown in inset); right panels: cooling 
sequence of the resistor when V{t)/V = 1 — H{t) (as shown 
in inset). Times are given in units of the electron-phonon time 
at temperature Te : r^-ph (Is) = 7/ (Sr|) . Blue dashed line 
is \/2 T* , red dashed-dotted line is 1 — 0.86 exp(— 4.2 r*), with 
r* = t/re_ph (Te) . 



(which corresponds to actual experimental conditions in 
Ref. 27,28 ). We now show that despite the short duty 
cycle d = tp/tr = 0.005, heating is not negligible. In 
our example, v ~ 30, therefore the relevant timescale 
when V is applied is ro_ph (Is) = 7/ (ST|) ~ 10 ^s. 
Equation (23) with 6lph = 15mK/Ts = 0.1 gives T{tp) = 
0.17 Ts = 25 mK, indicating slight heating by the first 
pulse. The resistance then cools down during a time 
tr before the next pulse is applied, following Eq.(21), 
to r = 0.169 Te, hardly less than at the end of the 
first pulse'^^. The temperature rises further during the 
next pulses, till steady oscillations establish. The full 
time evolution of T shown in Fig. 7 is obtained by iter- 
ating Eq.(21). At each pulse, the temperature rise gets 
smaller than during the preceeding pulse, because the 
starting temperature is larger and the heat transfer to 
phonons becomes more efficient. For the same reason, 
the cooling between the pulses gets more and more ef- 
ficient. At i > 250 /zs, a stationary regime is reached, 
with the reduced temperature oscillating between 6'inin 

to ^max such that /g 2"°'' dw/(H- - w^''^) = 2Tp and 

^^^2 d^f/(6'ph - w^''^) = 2(t; -t;). One obtains 6l„,i„ = 
0.3T(r„,i„ = 49 mK) and 0,nax = 0.36 (r„,ax = 54 mK). 
The amplitude of the oscillations A6' — 6'max — ^min is 
therefore very small. However, it increases with tp, as 
shown in the inset of Fig. 7, and can become sizeable. 

The main features of the time evolution of the temper- 
ature can be calculated more simply, from the average 
Joule power d x V'^ / R. Using section II, the character- 
istic temperature is then = 52 mK, which fits with 
the average temperature in the stationary regime of the 
pulse sequence. According to section IIIB, this tempera- 
ture is reached in a time To_ph(T£'^) ~ 220 ^s. The rise of 
temperature with time calculated with Eq. (21), in which 
all quantities (0, r*) are calculated using ^ is shown 
as a dotted line in Fig. 7, and reproduces well the overall 
behavior. The amplitude A6' of the temperature oscilla- 
tions can be evaluated using Eq. (23) under the assump- 
tion that the starting temperature is , which is a good 
approximation for oscillations of small amplitude, and 
considering that the voltage V is always present during 
the pulse of duration tp, i.e. with 6'(0) = /T^. = d^^^ 
and T* = tp/To_ph (Te) . If 6'ph «C 1 and c? ^ 1, one ob- 
tains A9 w 9{0)t*, an approximation only 20% larger 
than the exact result in the worst case of the inset of 
Fig.7 (tp^ 2 fis). 



IV. SUMMARY 

The solution of the heat equation in a resistor 
is determined by a characteristic temperature = 

(UV^r^i?)^^^ . If u = eU/fcsTE < 1, cooling by phonons 
is negligible and T{x) is given by Eq. (3), the average 
temperature by Tav = VS/S (eV/kB). if v > 10, the tem- 
perature is {T^+T^i^y^^ except at distances shorter than 
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^ ^Ljv from the ends. At w > 4 and Tph = 0, the av- 
erage temperature is Tav « — 1.16/w). Using these 
results, we have calculated the decay of the Fano factors 
Fi and F3 relative to the second and third cumulants of 
current fluctuations with the resistor length L. We have 
also addressed time-dependent situations to describe the 
heating and cooling of resistors. If < 1, the characteris- 
tic timescale is the diffusion time tjj and heating follows 



%/3 



eVy^l — exp(— lOr), cooling ksTg^vit) 



-^ey exp(— 5r), with r = t/ro- If v » 1, the instanta- 
neous relaxation time is Te_ph (T) =7/ (ET^) . Heating 
from T(0) to Ts is achieved in a time Tc_-ph , fol- 
lowing 0{t*) ', 



\ T* at short 



/ 0(0)2 + 2(1 + 61,^- em 

times and 0{t*) « 1 — 0.86 exp(— 4.2 r*) at long times, 
with 9 = T/Tt. and r* = t/'^o-ph (Te). Cooling from a 
temperature Tg occurs very slowly, along a powerlaw: at 

Tph = 0, r(i)/To = (1 + 3VTe-ph(To))"'/'. 

Finally, we recall that the actual temperature can be 
higher than the predictions made here for at least two 
reasons. First, the electronic temperature can be larger 
than Tph in the connecting wires because of their fi- 
nite resistivity-'^^'^^ or because of imperfect thermaliza- 
tion to the cryogenic unit. Second, we have neglected 
the Kapitza resistance^*^, due to which the phonon tem- 
perature inside the resistor can differ from the bath tem- 
perature Tph. However, this latest effect is relatively less 
important in very thin resistors because the ratio of the 
heat flow from electrons to resistor phonons to the heat 
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FIG. 7: (Color online) Main panel: Time-dependence of 
the temperature of a commercial macroscopic surface mount 
500 n resistor (see text) heated by voltage pulses (bottom 
curve) of length tp — 0.1 /^s applied every 20 fis. Dashed- 
dotted line is the predicted heating with the average Joule 
power. Inset: (blue dashed lines) minimal and maximal tem- 
perature reached in the stationary regime as a function of the 
pulse length tp. In red solid line, temperature evaluated with 
the average Joule power. 



flow from resistor phonons to substrate is proportional 
to the film thickness^. 



V. APPENDIX: CALCULATION OF THE 
THIRD CUMULANTS OF CURRENT IN 
PRESENCE OF ELECTRON-PHONON 
SCATTERING 

The calculation of the third cumulant of current in 
presence of electron-phonon scattering is an extension of 
the expressions of Pilgram et al.^^. The third cumulant at 
zero frequency is expressed as a function of the correlator 
between temperature and current fluctuations: 



^3 = ^ / 



(24) 



To calculate the integrand, we start from the stochastic 
diffusion equation for the fluctuations Sf of the electron 
energy distribution function 



(25) 



= - 1 A 

de L dx~ 

with 5Iee the linearized electron-electron diffusion inte- 
gral, 6Ie-ph the linearized electron-phonon diffusion in- 
tegral, SF^'^P and SF'^^ random extraneous sources asso- 
ciated with electron-impurity and electron-electron scat- 
tering. The correlation function of extraneous sources 
is 

{SF'"'P{e,x)SF'"'P{e',x')), ~ ^ 



. ^ S{x - x')5{e - e') 

X f{e,x)[l- f{e,x)]. (26) 

The energy distribution function is assumed to have 
a Fermi shape with coordinate dependent temperature 
T{x) and electrical potential (f>{x) : 



fie,x) 



1 + exp 



e(j){x) 



ksT^x) 



(27) 



To derive the correlator {5Te{x)6I)^, Eq. (25) is multi- 
plied by e and integrated over energy^^, assuming that 
the rate of energy dissipation associated with electron- 
phonon scattering is of the form 



(28) 



The electron-electron collision integral and the associated 
extraneous source drop out because of energy conserva- 
tion, and one obtains 



d_ _ 1 
dt Tpi dx'2 
1 9^ , 
Tp, dx^ ^ 



4 / dee-^6F"''P. (29) 
L ox 
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Now wc multiply Eq. (29) and the equation for the fluc- 
tuations of the total current, which in the low- frequency 
limit reads^^ 



61 



hi 



(30) 



Upon averaging, it gives in the low-frequency limit 

^ [LoT{5T{x)6I)a;] - LoaT^ {ST{x)6I)u. 



[m{x)dlU + 1^1 deef{l - /). (31) 



with a = bTiVlR/ Lq. The right-hand side of this equation 
was calculated in Ref. 17 . The solution of this equation 
may be written in a symbolic form as 



{5T{x)5I), 



2kf 



/ d{(t)T) 
1 dx 



dx"^ 



dx"^ 



dT 
dx 



(32) 



where the symbol {d^ /dx^ — /) ^ is the Green's function 
G{x,y) such that {d"^ jdx^ — f)G{x,y) = S{x — y) and 
G{0,y) = G{x,0) = G(l,y) = G(x, 1) = 0. Using = 
— Vx, the expression in brackets greatly simplifies; 



I dx 



dx'^ 



dx'^ 



dT 
dx 



y(T-2Tav). 



(33) 

To calculate the third cumulant of the current, one has to 
solve Eq. (32) and substitute the solution into Eq. (24). 
The generalized Fano factor = Ss/e'^I is then 



F3 



_ 36 /■ 

~^Jo ^T{x) \dx^ 
_ 36 /'V . 1 



' (^-aT\x) 



Ax Ay 
10 ^(x) 



{T - 2Tav} 
(34) 

Giie,x,y){e{y)-20,^} (35) 



with G\{6,x,y) the Green's function such that 



_ 15 

\ dx"^ TT^ 



2 /)3/ 



{x)]Gi{e,x,y)=5{x-y), (36) 



which can be calculated from^^ 
Go (a;, y) = [ ^ 1 = mm{x, y) (max(a;, y) - 1) (37) 



using 



G^ = {l-'^v^Goe^{x)] Go. (38) 



In practice, we performed this calculation by discretiza- 
tion of the coordinates and matrix inversion: the resistor 
is cut into A'' pieces of length s = 1/N, and the function 
Go (a;, y) is represented with a matrix G° such that 



G?. = -- min(i, j)(max(i, j) - 1) (39) 



(0 <i,j < N). The term i| GoO'-^ix) is represented by 
the matrix F build on the calculated temperature profile 
9{x) using 



(40) 



We then invert the matrix A with Aij = -5ij — Fij and 
compute G^ = A-\G°. Finally, 

36 ^GlA-0U:) + 20^„} 



e{i) 



(41) 
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